New Methods in Applied and Computational Mathematics (NEMACOM'98) 
Proceedings of the Centre for Mathematics and its Applications 



1 



Approximate Models of Dynamic Thermoviscoelasticity 
Describing Shape-Memory-Alloy Phase Transitions 

R.V.N. Melnik and A.J. Roberts 

Department of Mathematics and Computing, University of Southern Queensland, Toowoomba, 
QLD 4305, Australia. 

1. Introduction 

In this paper we consider two models for the approximate description of thermome- 
chanical behaviour of viscoelastic materials. Accounting for thermal fields in such a de- 
scription is important for all viscoelastic materials ranging from viscous fluids to elastic 
solids. The viscoelastic behaviour typically combines viscous and elastic properties and 
the relative proportion of this combination strongly depends on thermal characteristics 
of the material. Moreover, with changing thermal conditions, it is sometimes difficult to 
decide whether a particular material is a solid or a fluid. The key points in such decisions 
belong to the time of observation and to the choice of constitutive relations which couple 
stresses, deformation gradients, thermal fluxes and temperature. 

Our analysis is based on the nonlocal theory of continuum mechanics which considers 
constitutive variables defined at a point as a function of their values over the whole spatial 
domain of interest rather than as a function at that point only [2]. This approach of 
rational mechanics allows us to derive a general model that is suitable for the description 
of thermomechanical behaviour of materials under a wide range of temperature and 
loading patterns. In our models we allow for the dependency of stresses not only on 
the deformation gradient and temperature but also on the rates of their changes. Such 
considerations put us closer to real situations where the time-dependent coupling between 
temperature and stresses have to include the velocity of the deformation gradient and 
the speed of thermal propagation. Another novelty of our paper is the accounting for 
finite speeds of thermal disturbances. We define the constitutive relationship for thermal 
fluxes using the Cattaneo-Vernotte equation which includes the classical Fourier law as 
a special limiting case (in the limit of zero relaxation time for heat fluxes). In particular 
this approach is critical in modelling short transient states in low temperature regimes. 

During recent years a number of papers were devoted to the development of mathemat- 
ical theory of thermomechanical phase transitions (see [19, 27, 12, 13, 1] and references 
therein). The majority of those papers dealt with important theoretical issues of models 
such as well-posedness and the global asymptotic behaviour of solutions. However, only 
a few papers have been devoted to the description of computational results using those 
models (see, for example, [20, 15] and references therein). Almost all developed models 
take into account neither the rate of thermal disturbances nor the relaxation time of 
thermal fluxes. However, the importance of these issues are well known in dynamic hy- 
perbolic thermoelasticity where mathematical procedures and computational techniques 
have a longer history compared to that in thermoviscoelasticity [16, 23]. 

In dealing with the three main physical quantities of continuum mechanics (stresses, 
deformation gradients and displacements) it is important to take into account their cou- 
pling to the thermal field. This allows us to construct efficient mathematical models 
for the description of complicated phenomena, such as hysteresis, which are becoming 
increasingly important in a wide range of applications. In this paper we apply the de- 
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veloped models to the description of shape memory alloy effects in a large bar. It is 
well-known that for many types of shape memory materials the dependency of stresses 
on the deformation gradient upon loading and unloading is significantly different. Ap- 
plying a large load at a low temperature, we may get a residual deformation gradient, 
which typically vanishes upon heating. The restoring of the original shape is referred to 
as the shape memory effect. This effect is discussed with two numerical examples. 
The rest of the paper is organized as follows. 

• Section 2 provides the reader with basic preliminaries and notation. 

• The general formulation of the model is given in Section 3. In this section we specify 
the model for internal energy and derive restrictions on the model imposed by the 
second law of thermodynamics. 

• In Section 4 we incorporate the Cattaneo-Vernotte equation for heat conduction 
into our model. 

• Section 5 deals with the Landau-Devonshire model for the free energy function. 
The constitutive relation connecting stresses and the deformation gradient is also 
discussed in this section. 

• In Section 6 we consider a one-dimensional model of thermoviscoelasticity and dis- 
cuss the consequences of non-convexity of free energy function. 

• Some numerical results are presented and discussed in Section 7. 

• In Section 8 we use centre manifold theory to derive an approximate mathematical 
model for the description of thermomechanical behaviour of viscoelastic materials. 

2. Preliminaries and Notation 

Assume that an object of interest (a solid, fluid, gas or plasma) occupies the volume 
V in a fixed reference spatial configuration Q at a certain time to- This object in its 
spatio-temporal configuration will be referred to by the generic name "system". We 
aim to develop an efficient mathematical description of the dynamic thermomechanical 
behaviour of the system. 

Let x = (xi,X2,xs) be material (Lagrangian) coordinates of a material point of the 
system in the configuration Q at time t . Then the dynamics of the system is determined 
by the spatial displacements u = (ui, 112,113) of such material points as a function of the 
reference position, x, and the time of interest, t. The partial derivative of displacement 
with respect to x is identified with the symmetric strain tensor 



and the time derivative of the function u is identified with the velocity of the system 



e = sym 



<9u(x, t) 
<9x 



or €ij (x, t) 



2 dxj dxi 



1 dv,i(x.,t) duj(x,t) 



i,j = 1,2,3, (2.1) 




du 



or 




(2.2) 



In (2.1) we require that det(J + e) > which precludes a possibility of compression of the 
matter to zero and guarantees the local invertibility of x + u(x, t) [22]. Since the time 
derivatives are understood in the Lagrangian sense, x is kept fixed in (2.2). 
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3. Nonlocal Models of Thermoviscoelasticity 

The equation of motion requires information on forces acting per unit area of the 
matter and, hence, in a natural way, involves the concept of stresses. The stress is not a 
mere function of of the deformation gradient, as it is often assumed. It also depends on 
temperature of the matter, its rate of change in time and the rate of change of deformation 
gradient e. Let po(x, to) > be the density of the matter (the mass per unit volume) in 
the reference configuration Q at time to and p(x, t) be the density of the matter at time 
t where t — to is sufficiently small. Then, in the Lagrangian system of coordinates (x, t), 
the equation for balance of mass is written in the form [22] 

p(x,t)det(/ + e(x,t)) =po(x,t ). (3.1) 
The equation of motion has the following form 

p^ = V x -s + F with F = p(f + f)-pv, (3.2) 

where f is a given body force per unit mass, p and f are nonlocal mass and force residuals 
respectively, and s is the stress tensor. 

In Lagrangian coordinates the equation for energy balance has the form 

^( e + y) -V x -(s-v) + V-q = p^ + ^ + f - v "^( e + y))> ( 3 - 3 ) 

where e is the specific internal energy of the system, v 2 = v • v, h is the heat source 
density, h is the nonlocal energy residual (see [2] for conditions on localised residuals) 
and q is the heat flux. The scalar multiplication of (3.2) by v gives 

p-^- — v • (V • s) = (F, v) = p(f + f ) ■ v — pv 2 . (3.4) 
Taking into account normalisation, from (3.3) and (3.4) we get 

p^-s T :(Vv) + V-q = s, (3.5) 

3 

where a : b = a^bij is the standard notation for the rank 2 tensors a and b and 

g = p(h + h)-pf-v-p(e-?Pj. (3.6) 

The right-hand sides of equations (3.2) and (3.5) incorporate into the model nonlocal 
and dissipative effects of thermomechanical waves. As we shall see in the next section, 
under appropriate constitutive relations it is also possible to allow for a relaxation time 
for acceleration of the motion in response to applied gradients such as the deformation 
gradient and the temperature gradient. 

We assume that there exists a one-to-one entropy function of the system state. We 
denote the density of such a function by rj, and then the second law of thermodynamics 
is 
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where £ is the entropy source density, r is the entropy flux density and £ is the nonlocal 
entropy residual. 

The system of equations (3.2), (3.5) combined with inequality (3.7) provides the gen- 
eral mathematical model for the description of thermomechanical behaviour of dynamic 
systems. The macroscopic modelling of such systems starts from the choice of constitu- 
tive relationships. We assume the existence of a functional ^ invariant under a time shift 
and chose this functional in the form of the Helmholtz free energy 

V = e-9ri, (3.8) 

where 9 is the temperature of the system {9 > 0, inf( Xjt ) 9 = 0). We also assume specific 
forms for the entropy flux and the entropy source density as 

r = q/9, S = h/0. (3.9) 

Using (3.8) in (3.5) and taking into account that 

V-q = 0V-(q/0) + (q-V0)/0, (3.10) 

from (3.7) and (3.9) we get the nonlocal formulation of the Clausius-Duhem inequality 

+ - f .v-^-(^„o. ( ,n, 

The latter inequality together with requirements on localisation residuals (see [2] for 
details) impose restrictions on the choice of nonlocal residuals and the functions rj, s and 
q. We assume that the entropy density is given in the form 

"=-»■ (3 ' 12) 

Finally, we have to specify the constitutive relationships that couple stresses, deformation 
gradients, temperature and heat fluxes 

$i(s,e) = 0, $ 2 (q,0) = O, (3.13) 

where it is implicitly assumed that these relations may involve spatial and temporal 
derivatives of the functions. In Section 4 and 5 we specify particular forms for $i and 

$2- 

4. The Cattaneo-Vernotte Model for Heat Conduction 

The choice of the function $ 2 in (3.13) is made using the Cattaneo-Vernotte model 

ci + r ^ = -k(9,e)V9, (4.1) 

where To is the dimensionless thermal relaxation time and k(9, e) is the thermal conductiv- 
ity of the material (typically k = 1 + /39 with the given dimensionless coefficient f3). Such 
a choice is made in order to account for the finite speeds of thermal wave propagation and 
thermally induced stress wave propagation coupled to the deformation gradient [11, 18]. 
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In order to incorporate equation (4.1) into the general model of thermoviscoelasticity we 
use a consequence of (3.5) 

^_.| [sr:( Vv) ]+ .V.(|)=.|. (4,) 

On the other hand, from (4.1) we get 

V • q + t V • (^j = -V • (W6) . (4.3) 

Then from (3.5), (4.2), (4.3) we obtain the energy balance equation in the form 

p m + pr °w ~ sT : (Vv) " r °^ [sT : (Vv)] " v ' {kvd) = G ' (44) 

where 

G = g + r ^. (4.5) 

During recent years, the interest in such a hyperbolic approach in the analysis of materials 
with memory has increased [6]. 

5. The Landau-Devonshire Model for the Helmholtz Free Energy and the 
Stress-Strain Relation 

We start from the consideration of the one-dimensional case assuming the following 
approximation for the free energy of the system 

(0, e) = Vo(0) + M0)Me) + (5-1) 

where ipo{0) models thermal field contributions, ipi{9)ip2{^) models shape-memory con- 
tributions and ip3(e) models mechanical field contributions. These models are chosen in 
the following forms 

ip (9) = a - axOlnO, ipi{6) = \a 2 6 , tp 2 {e) = e 2 , 

1 11 ( 5 - 2 ) 

M e ) = -^ 2 0ie 2 - -a^e 1 + -a 6 e 6 , 

where all ctj and 6\ are positive constants. The model (5.1)-(5.2), known as the Landau- 
Devonshire model for the Helmholtz free energy, covers a number of important practical 
cases. However, it belongs to the class of models which is difficult to investigate compared 
to the Landau-Devonshire-Ginzburg model. In the latter case an additional "smoothing" 
term in (5.1), known as the Ginsburg term ^u xxxx ^ allows us to obtain a bound of the 
deformation gradient (strain) using a well established technique [5] . 

Remark 5..1 A number of important characteristics of phase transformations (such as 
the size of hysteresis) may depend on the contributions of the interfacial energies. These 
contributions are often modeled with the Ginsburg correction term. However, the Ginsburg 
coefficient can only be determined in approximate order [28] and in the general case 
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this coefficient may not be temperature-independent. Another way to account for the 
contributions of inter facial energies is to to take the free energy in the form [4, 17] 



* = (1 - z)&(e, 9) + # 2 (e, 9) + z(l - z)4> 3 , 



(5.3) 



where z is the volume fraction of martensite (i.e. the product phase), (1 — z) is the volume 
fraction of austenite (i.e. the parent phase), ipi, ip2 ^re the free energies of austenite and 
martensite respectively and is the contribution from the interaction effect between 
austenite and martensite. We will not pursue these ideas in this paper. 

Remark 5. .2 Some authors include a linear term a°9 into ipo(0). This term has no 
bearing on the final model and changes only the value of the coefficient of 9 in the internal 
energy representation (see formula (6.2)), and thus is omitted. 

In the general case for the choice of the function $i in (3.13) we allow the dependency 
of the stress on the rate of temperature and the deformation gradient 

'89 de ~ 



p(9,e) + \ 



dt ' dt 



where 



tn N df A (89 de\ _ tn .de .89 



(5.4) 



(5.5) 



It is straightforward to deduce 
p(9, e) = a 2 9e + 



<9^ 3 (e) 
de 



a 2 e{9 — 9i) — a 4 e 3 + a 6 e £ 



(5.6) 



One-Dimensional Hyperbolic Approximation of Shape-Memory-Alloy 
Dynamics 

Using the model (5.1) and (5.2), from (3.12) we get 



rj = Q!i(l + ln^) — -ct^e 2 



(6.1) 



This enables us to find the internal energy of the system as a sum of thermal and me- 
chanical fields contributions 



ciq + ct\9 



1 1 1 

-a 2 9\e a 4 e H — a 6 e 

2 4 6 



a + ol x 9 + ^(e) 



(6.2) 



The substitution of (6.2) into (4.4) leads to the final form of the energy balance equation. 
In particular, assuming symmetry of the deformation gradient tensor, we get 



pot\ 

where the meaning of A is 



09 d 2 9 



+ A(e, 6) - V • (kV9) = G , 



M n\ in 9e 8 

A i e , °) = -P«2 < 9e— + r — 



(6.3) 



d_ 

dt 



'de 
dt 



d9 \ . ,de d 



dt 



dt 



dt 



de 
dt 



(6.4) 
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Equation (6.3) is solved together with the equation of motion (3.2) with respect to (u, 0): 



dO 



d 



( d 2 U \ 



dt 2 

2 



dt \ dtdx 



du d 2 u 
dx dtdx 

d6 d 2 u 
dt dtdx r ° dt 



d ( 'du d 2 u \ 



dt \ dx dtdx 



d 



d 2 u 



d_ 

dx 






(6.5) 



- p 



— v- 



dx 2 dt dxdt 



where C v = pa±, k\ = pa 2 , k 2 = pa^, k 3 = pa 6 , p = pp, v = pu. 
The initial conditions for the model (6.5) are chosen in the form 



u(x, 0) 



u°(x) 



u 



(x); 0(x,0) = 0°(x) 



e\x) 



(6.6) 



for given functions u°, u 1 , 0°, 1 . There are several distinct choices for boundary condi- 
tions to be used in our computational experiments. Mechanical boundary conditions are 
taken in one of the following forms (L is the length of the structure): 

• "stress-free" boundary conditions: s(0, t) = s(L,t) = 0; 

• "pinned end" boundary conditions: u(0, t) = u(L,t) = 0; 

• or mixed mechanical boundary condition: s(0, t) = 0, u(L,t) = 0. 

When displacements are given on boundaries, a priori bounds on strains are generally 
unknown which complicates the mathematical analysis of the problem. Computational 
results presented in Section 7 deal with this case. Thermal boundary conditions are 
chosen in one of the following form 

• "thermal insulation" boundary conditions: q(0,t) = q(L,t) = 0, which reduce to 
d6 d6 

_(0,t) = ~-(-M) = for the Fourier law (see (8.2)); 

dO dO 

• "controlled flux" boundary conditions: 7p(M) = 0, -k—(L,t) = (3[0 - 0°(t)]; 

• or fixed temperature ("uncontrolled energy flow") boundary conditions: 0(0, t) = 
0(L,t)=0. 

In the last case additional assumptions are needed. By using the Leray-Schauder principle 
we have analysed the Cauchy problem for nonlinear hyperbolic model of thermoviscoelas- 
ticity (6.5). Our procedure makes use of the Lumer-Phillips theorem and the technique 
developed in [5]. We shall address details of this technique elsewhere. 

Our final remark in this section goes to the choice of the function \1/ in the form (5.1) 
and (5.2) that brings major difficulties in the investigation of the model (6.5). Strictly 
speaking, the free energy function strongly depends upon the statistics of the phenomenon 
and has to be derived from a statistical model. Since van der Waals work on statistical 
mechanics it is a common practice to choose this function as a non-convex function of 
e [14]. When dealing with shape memory alloys, minima of this function are known to 
correspond different phases of the material. For example, in the case of three minima, 
we expect one austenitic and two martensitic phase (see, for example, [8, 19, 28, 12]). 
Temperature plays a crucial role in the phase transition. Depending on the value of 
temperature, the material may alternate between a single thermodynamically unstable 
nonmonotone branch and multiple unstable branches. The character of this instability 
depends not only on the deformation gradient and temperature, but also on the rates of 
their changes. 
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7. Computational Experiments 

In this section we present some numerical results on the thermal and mechanical control 
of a rod (L = 1cm) with a shape-memory-alloy core. The parameters of the Cu-based 
core are taken as follows 

k = 1.9 x l(r 2 cmg/(ms 3 K), p = ll.lg/cm 3 , C v = 29g/(ms 2 cmK), 9 l = 208K , 
fci = 480g/(ms 2 cmK), k 2 = 6 x 10 6 g/( ms2cmK )> h = 4.5 x 10 8 g/(ms 2 cmK) . 

We use model (6.5) with r = = /! = z/ = 0, initial conditions (6.6) and "pinned end 
& controlled flux" boundary conditions. This model was straightforwardly reduced to a 
differential-algebraic system in tz = (u, v, 9) T and stress s using second-order accurate 
spatial differences on staggered grids: 

D _ = r, . = *,(«_«,)__*, +* 3 (-j , (7.i) 

where D is the diagonal matrix with diag(D) = (1, p, C v ), f = (/i, /2, /s) T and 

/i = «, A = l £ + ^, /s = ^ + M^ + G. (7.2) 
OX ox 1 ox ox 

The developed code is robust and much simpler compared to computational procedures 
previously reported in the literature for shape-memory alloys [20, 15, 13]. 

Experiment 1 (thermal control of phase transformations). In this experiment 
we set uniform forcing F = 500g/(ms 2 cm 2 ) and vary heating conditions given by G = 
3 757rsin 3 (t7r/6)g/(ms 3 cm). We assume that the initial displacements are given in the 
form 



u°(x) 



-0.11809a;, < x < 1/6, 

0.11809(x-l/3), 1/6 < x < 1/2, , . 

0.11809(2/3 -x), l/2<x<5/6 ^ ' 

0.11809(ar- 1), 5/6 < x < 1 



and take the initial temperature as 9° = 200K. Figure 1 (obtained with time step 
7 x 10~ 4 ms and space step l/24cm) demonstrates the transformation of 2M + + 2M~ 
martensites into an austenite (visible in the region of zero strain and displacements with 
superposed elastic vibrations as seen most clearly in the velocity field) after sufficient tem- 
perature has reached. Then upon cooling we observe a first order (martensitic) transition 
from the high temperature phase (austenite) to the low temperature phase (martensite). 
Upon the return to the low temperature regime the stable attractor with this applied 
thermomechanical forcing is not the original configuration but only two distinct marten- 
site phases. The transformation [2M + + 2M~] — > A is accompanied by a decrease in 
temperature whereas the transformation A — * [M + + M~] is accompanied by an increase 
in temperature. 

Experiment 2 (mechanical control of phase transformations. In this experi- 
ment we set G — 0, but vary the mechanical loading according to F — 7000 sin 3 {id/ 2) 
g/(cm 2 ms 2 ). Starting from the austenite configuration {u° = 0) at intermediate temper- 
ature 9° = 255K we observe (see Figure 2 where the time step was 8 x 10~ 4 ms and the 
space step was 1/1 6cm) the transformation A -> [M+ + M~] —> A —> [M~ + M+] -> A. 
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FIGURE 1. Thermally induced phase transformations. 



In this experiment we observe the almost immediate transformations of austenite phases 
into two martensites upon the increase/decrease in loading. Note the relatively large 
heating/cooling associated with the transition into/out of martensite phase. A similar 
behaviour under different thermomechanical conditions was also observed in [20, 15]. In 
our code we have also incorporated the Ginsburg term by adding r yu xxxx to f 2 in (7.2). 
With reported values of the Ginsburg coefficient (7 ~ 1CT 10 — 1CT 12 ) the Ginsburg term 
has a negligible effect on the thermomechanical behaviour of shape-memory alloys in the 
group of experiments described here. Accounting for interfacial energy contributions and 
the influence of mechanical and thermal dissipations on the dynamics of memory material 
require further investigation. 

8. Construction of Approximate Models for Dynamic 
Thermoviscoelasticity Using Centre Manifold Theory 

The model described in Section 6 will provide a good approximation of thermomechan- 
ical behaviour of a large shape memory alloy bar (see applications in [3]) only in the case 
the bar can be modelled by a thin rod with a shape memory alloy core. As an alternative 
to that model, in this section we construct a new model which is derived directly from 
the 3D model for shape memory alloy evolution (see (3.2), (6.3)) using centre manifold 
techniques (see, for example, [24]). 

We assume that the shear stress in equation (3.2) is determined by its three compo- 
nents, the quasi- conservative component, s q , the stress component due to mechanical 
dissipation, s m , and the stress component due to thermal dissipations, s*, (the latter is 
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Figure 2. Mechanically induced phase transitions. 



assumed to be negligible at this stage) 

s = s 9 + s m + s', with s q = p^-, s m = p/i^, s' = 0. (8.1) 

In the general case the heat flux is determined as the solution of equation (4.1). An 
approximation to this solution is provided by the following generalised form (see [19] and 
references therein) 

q = -kV6-a^-, a>0, (8.2) 

which we will use with a = when (8.2) turns into the classical Fourier law. 
The internal energy function e is defined from (3.8) by 

e = *-» w (8.3) 

In order to complete the formulation of the problem we specify a model for the free 
energy function However, in the general 3D case one cannot use the shear strain as 
the order parameter as we usually do for the ID case. One of the first approaches to deal 
with the 3D challenge was the Fremond model. This model uses different expressions for 
free energy functions for different phases (see, for example, [10]). All these expressions 
are essentially of the Landau-Ginsburg-type and contain the term 7/2Vtr(e) with 7 > 
introduced in order to smooth possibly very sharp spatial phase separation. In this 
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paper we use a different approach proposed in [9]. This approach generalises the classical 
Landau-Devonshire-Falk theory for shape memory alloys to the 3D case. The free energy 
function, based on the expansion up to sixth order in a single shear strain component [8] , 
was extended to the three-dimensional case using the group theoretical approach proposed 
in [21]. In contrast to some other models (see, for example, [10]) strain-gradient terms 
are not involved in his expansion. We make use of this expansion and apply the following 
general representation of the free energy function 



(8.4) 



where independent material parameters of the ra-th order for n 
through strain invariants, X™ , as follows 



1, 2, ... are determined 



.5) 



The upper limit of the sum in (8.5), j n , is the number of all invariant directions associated 
with a representation of the 48th order cubic symmetry group of the parent phase (see 
details in [9]). In order to adequately describe thermomechanical behaviour of shape- 
memory alloys we need to account for 6 terms in the sum of the expansion (8.4). In this 
case we have to determine 32 material parameters that make the application of formulae 
(8.4)-(8.5) fairly complicated. Using physically justified assumptions it is possible to re- 
duce the number of required parameters. To achieve this, we make the same assumptions 
as in [9]. They conclude that odd degree invariants can be neglected in the expansion. 
Taking invariants up to the sixth order results in a representation with only 10 material 
constants which may depend on temperature 



# = ip°(o) + j2 m + E ^ + E $3 

j=i j=i j=i 



(8.6) 



(we do not neglect the contribution of ijj°(9) as was done in [9]). The strain invariants Xf 
of second, forth and sixth orders of the 48th order cubic symmetry group of the parent 
phase are 

A = ^(tr(e^)) 2 , X 2 = ^(2e 33 - e n - e 22 ) 2 + ^(e n - e 22 ) 2 , 



X 2 



-23 



+ e z 13 + e 



12) 



Z\ = (X 2 ) 2 , 



X 4 

x 2 



: 23 



+ e? 3 + e' 



12) 



It = (X 2 ) 3 



7"4 _ /-7-2\2 7-4 _ 7-27-2 7-. 
- L ?> ~ V X 3/ > X 4 — X 2 X 3> X f 



6 23 



77(2633 — en — e 22 ) — -(en — e 22 ) 

D 2 



+ 



: 13 



i(2 £33 



en 



e 22 ) + -( e n 



£22 



+ g e i2( 2e 



33 



e ll — e 22) 2 ) 



(8.7) 



X 6 

x 2 



7777(2633 — e n — e 22 ) 2 ( — (2t33 — e±i — e 22 ) 



36 



36 



4^11 



e 22 )' 



The ten material constants ?/>™ in (8.6) differ from alloy to alloy and we use coefficients 
derived for Cu-based alloys [9] (units used here are consistent with those used in Section 
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7 for our numerical results on the dynamics of shape- memory alloys): 

V>i = 5.92 x 10 6 g/(ms 2 cm), ^2 = ( L41 x !0 5 + 46(0 ~ 300)) g/(ms 2 cm), 

^3 = (1.48 x 10 6 - 940(0 - 300)) g/(ms 2 cm), ^° = -ai0rn[(0 - 9 )/9 ] g/(ms 2 cm), 

^ = (_i.i82 x 10 8 + 3.55 x 10 5 (9 - 300)) g/(ms 2 cm), 

V>2 = 3.13 x 10 9 g/(ms 2 cm), ^ = 1.64 x 10 9 g/(ms 2 cm), (8.8) 

$\ = -5.53 x 10 s g/(ms 2 cm), ^ = -4.27 x 10 8 g/(ms 2 cm), 

V>i = 3.35 x 10 10 g/(ms 2 cm), ^2 = 3.71 x 10 11 g/(ms 2 cm). 

Other material parameters are taken to be the same as those given in Section 7. We 
are interested in the construction of an adequate model for the description of thermo- 
mechanical behaviour of thin slabs in shape memory alloy materials. Starting from the 
3D Falk-Konopka model and using centre manifold techniques (see, for example, [24]) we 
derive systematically an accurate low-dimensional model for the dynamics of the slab. 
The shape memory alloy is assumed to be of very large extent in the x = x\ direction 
compared to its thickness of 26 in the y = x<i direction (—b < y < b). For the sake of 
convenience we use a new temperature variable & = 9 — 9$ where here 9q = 300. For 
simplicity of the analysis we assume zero dissipation, a — /i — 0, and that there is no 
motion nor dependence in the x 3 direction. 

A model for the dynamics of modes which vary slowly along the slab is derived for 
the unforced dynamics, F = 0, G = 0, and when "zero-stress & thermal-insulation" 
boundary conditions are specified on y = ±b. The derivation of boundary conditions in 
the "long" direction x requires a quite delicate analysis [25] and these issues will not be 
address here. We only note that "pinned & insulating ends" boundary conditions may 
be used as a leading approximation. Modelling the long-wavelength, small-wavenumber 
modes along the slab, we neglect all longitudinal variations and look for eigenvalues 
of the cross-slab modes. It can be shown that generally there is a zero eigenvalue of 
multiplicity five and all the rest are pure imaginary (as dissipation has been omitted). 
Thus there exists a sub-centre manifold based upon these five modes (see [26] for an 
existence theorem), called a slow manifold as these five modes evolve slowly. Note that 
being on a sub-centre manifold the models we construct only have a weak assurance of 
asymptotic completeness, see the discussion in [7]. The zero eigenvalue of multiplicity five 
corresponds to longitudinal waves, large-scale bending, and one heat mode. The leading 
order structure of the critical eigenmodes are constant across the slab. Thus letting an 
overbar denote the y average, the amplitudes of the critical modes are chosen in the form 

Ui(x,t)=m, Vi(x,t)=Vi, &(x,t)=W. (8.9) 

The low- dimensional model below is written in terms of these parameters. 

The construction of the low- dimensional model is based upon the ansatz that there 
exists a low- dimensional invariant manifold upon which the amplitudes evolve slowly: 

Ui = u i (u,v,e'), Vi = Vi(u,v,e'), e = r(u,v,e'), (8.10) 

BU 8V <9B' 

where — 1 = V, , — 1 = 9i (U, V, 9') , — = g e (U, V , 6') . (8.11) 
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This ansatz is substituted into the differential-algebraic equations of 3D thermo-visco- 
elasticity (3.2), (3.5) and is solved to some order in the small parameters d x , E = \\U X \\ + 
||V X || and •& = ||0' 1 1 with the computer algebra package Reduce. Thus, here we treat the 
strains as small, as measured by E, while permitting asymptotically large displacements 
and velocities. The displacement and temperature fields of the slow manifold, in terms 
of the amplitudes and the scaled transverse coordinate Y = y/b, are approximately 

ui « E/i -YbU 2x + 0.15(3 Y 2 - l)b 2 U lxx , (8.12) 
u 2 « U 2 - (0.9 - 3.05e-5 Q')YbU lx + 0.15(3 Y 2 - l)b 2 U 2xx 

- 141 YbU lx 3 + 1.00e-4 (3Y - Y 3 )b 3 V lx 2 U lx , (8.13) 
9 « 300 + 0' - 2.43e6 (3Y - Y 3 )b 3 (V lx U 2xx + U lx V 2xx ) 

- 25.1(7- 30 Y 2 + 15 Y 4 )V lx 3 U lx . (8.14) 



These expressions have errors O [E 5 + d b J 2 + ^ 5 / 2 J where the notation O (E p + d q x + $ r ) 

is used to denote terms involving d\E a {} c such that a/p + b/q + c/r > 1. The mechanical 
and thermal field approximations represented by (8.12)-(8.14) have cross-slab structure. 
In particular, the sideways deformation u 2 (which is a nonlinear function of the longitu- 
dinal strains) of the shape memory alloy feed back at higher order to contribute to and 
complicate the longitudinal and thermal dynamics. 

The model for the longitudinal dynamics on this slow manifold is 

dV 

p-^- = 2.97e6U lxx + 8.03e5 b 2 U lxxxx 

+ d x [(922 0'-O.O145e /2 )f/i :c - (4.28e9 - 1.31e7 Q')U lx 3 + 7.12ell U lx 5 

+ (2820 - 8.80e')6V 1:r 2 C/ 1:r +l.2Ab A V lx A U lx - 5A2e4b 2 V lx 2 U lx 3 

+ 0(E 8 + d 4 x +$ 4 ) . (8.15) 

The first line in the right-hand side of (8.15) describes linear dispersive elastic waves along 
the slab, whereas the second line gives the temperature dependent quintic stress-strain 
relation of the shape memory alloy. Since V lx = Ui xt , the remaining lines show effects 
upon this stress-strain relation due to rates of change of the strain. 

Note that to this order of truncation there is no coupling to the bending modes of the 
slab which to the same error is simply the beam equation 

= -9.91e5b 2 U 2xxxx + 0(E s + d 4 x + $ 4 ) . (8.16) 

There exists nonlinear coupling between the modes at higher order. 
The corresponding energy equation for the temperature is 



~dt 



C v ^ = kB' xx + (2.77e5 + 914 0' - 9.25 & 2 )U lx V lx 



+ (3.94e9 + 1.26e7 Q')V lx U lx 3 - (57.3 + 0.0117 Q')b 2 V lx 3 U lx 

+ 1.68el2V lx U lx 5 - 1.58e6b 2 V lx 3 U lx 3 - 0.0203 b 4 V lx 5 U lx 

+ lMe4b 2 U lxx V lxx + 9.22e4b 2 U 2xx V 2xx + d 2 x [-8151 b 2 U lx V lx 

+ 0(E 8 + d 4 x + # 4 ) . (8.17) 
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The first line in (8.17) describes the diffusion of heat generated or absorbed by mechanical 
strains, ®U lx Vi x . However, in the thin slab the internal pattern of strains causes a much 
more complicated distribution of heating and cooling as summarised by the remaining 
lines. It is expected that virtually all of these should be retained in order to be con- 
sistent with the quintic stress-strain of the longitudinal wave equation. Computational 
experiments with the model derived in this section will be presented elsewhere. 
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